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1. Introduction 


Since the paper of Nielsen and OlesenjlJ describing the vortex solutions of the Abelian 
Higgs model, many variations of this structure have been studied, giving rise to diverse phe¬ 
nomena. Superconducting cosmic strings |3] are a dramatic example of this. Cosmic strings 
arise as classical solutions of many Grand Unified Theories (GUTs) and other models of 
fundamental interactions. They may play an interesting role in cosmology since they are 
formed in the early universe, when the order parameter of spontaneous symmetry breaking 
is uncorrelated at large distances; later, they may perturb the space-time metric |§ [£f| and 
contribute to the density fluctuations influencing galaxy formation. In most previously 
studied variations of the basic Nielsen-Olesen vortex solution the cylindrical symmetry is 
preserved. We explore another variation in the basic structure of cosmic strings arising 
in some non-minimal extensions of the standard model, in which additional fields appear 
in the classical equations and the resulting configurations break the azimuthal symmetry. 
(These additional fields may arise, for instance, in an extended Higgs sector, or may be 
part of the matter sector of the theory.) We will discuss two classes of models: one in 
which pairs of string are confined due to a spontaneously-broken discrete gauge symmetry, 
and another class in which independently stable cosmic strings form bound pairs of nearly 
parallel strings because of short-range attractive forces. 

Although there are several mechanisms for obtaining strings with an asymmetric core, 
the macroscopic classical motion of the strings can be described by a geometrical effective 
action, and only the relative size of the couplings depends on the details of the short- 
distance structure. For this reason, we discuss here the conditions for obtaining binary 
strings, and report on the effect of this additional structure on the macroscopic behavior 
of the cosmic strings in a separate paper |Q. 

A simple physical argument will serve to illustrate the necessary condition for a bound 
state to occur. Recall [[6] that in the Nielsen and Olesen solution, two vortices (that is, 
cross-sections of a cosmic string) will attract or repel each other depending on the relative 
strength of the electric charge and the scalar self-coupling. For /3 = ()yy) 2 < 1, two n — 1 
vortices have a short range attraction, so that the ground state in the two-vortex sector 
is the n — 2 vortex. For j3 > 1, the magnetic force dominates, leading to a repulsive 
short range force. (The critical [|6| case f3 — 1 has additional symmetries, but it will not 
concern us here.) In the case of GUT cosmic strings, we will show that for some ranges of 
couplings the force between vortices is attractive at long distances but at short distances 
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the force becomes repulsive. This can easily happen, for instance, when the longer-range 
force is confining, that is, when the discrete symmetry (the holonomy of a single string) is 


spontaneously broken. This leads to a long-range force confining pairs of vortices; if the 
short-range interaction is repulsive (/? > 1 ), they do not coalesce into a single rotationally 
symmetric vortex. Such a structure can also arise when the balance between short-range 
repulsive forces and short-range attractive forces between two vortices has a local minimum, 
leading to a finite binding energy per unit length. In either case the classical ground state 
consists of a pair of finitely-separated vortices, or in the (3+l)-dimensional case, a bound 
pair of cosmic strings. We will then discuss examples of both realizations of binary cosmic 
strings, beginning with the confining case. 

2. Confined binary strings 

We will present two models in which binary strings arise because of long-range con¬ 
fining forces. In both cases we have ‘unary’ cosmic strings formed at one scale, which are 
charged under a discrete gauge symmetry. At a lower scale, another scalar field breaks this 
discrete symmetry, and the original strings become boundaries of domain walls. These 
domain walls serve to bind together pairs of the original strings into binary strings, which 
are then in a sense quasi-one-dimensional domain walls, whose width and thickness are 
comparable in size. Later in the paper we will describe some GUT realizations of such 
models, such as the domain-wall-connected strings previously discussed in the literature 


i. 


In the first model, we consider a spontaneously broken U( 1) gauge theory with two 
charged scalars (ft and y; the field eft has twice the charge of y. The Lagrangian is: 


L = + \D 1 <p\ 2 + \D 2 x\ 2 - V. 


( 2 . 1 ) 


where D lfl (ft = (<9 M - ieA^)(ft, D 2 ^x = - |ieA M )y and 



We also define 



( 2 . 2 ) 


We shall show that for a range of values of 7 and the coupling constants Ai and n. and 
for all values of A 2 , the lowest-mass stable configurations are binary strings. The finiteness 
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of the mass per unit length of local strings (in contrast to that of global strings) derives 
from the proper long-range behavior of the gauge field; to ensure this, the winding number 
v 0 of the phase of 0 must equal twice the winding number of the phase of : iy 0 — 2v x . 
The smallest finite-mass solution then has — 2. In the absence of y there would be 
stable v<j, = 1 strings, but in this model separating the = 2 string into two of the 
original strings with = 1 creates a domain wall, and for large separations requires a 
linearly-increasing energy. We then need to determine the form of the — 2 string. 
We will demonstrate that for the range of parameters mentioned above, the cylindrically 
symmetric configuration is an unstable solution of the classical field equations; since we 
have argued that the string cannot separate, this indicates that the binary string is the 
lowest-energy classical solution. 

We will first demonstrate the instability of > 1 strings (and hence existence of 
binary strings) when we set the coupling k — 0 , and later demonstrate that this restriction 
is inessential. To start, consider the ((f), A) system if we totally ignore (i.e., decouple) 
the field y; then f3 — ( 777 ) 2 = (^r), and for convenience we define a — (3 — 1. In that 
case the proof proceeds along the lines of ||, where it was shown that cylindrically sym¬ 
metric n — 2 strings in the theory with a > 0 (and only one Higgs field, i.e. here we 
ignore y), are unstable against specific perturbations which correspond to the splitting of 
the string. The instability occurs because for a > 0 there are asymmetric perturbations 
around the cylindrically symmetric solution that lower the energy; for a — 0 these asym¬ 
metric perturbations of the cylindrically symmetric classical solutions are zero modes of 
the Hamiltonian. 

If we make the same perturbations and include the effects of the lighter Higgs field y, 
(so A 2 and 7 are nonzero), there will be additional positive contributions to the energy, and 
the line of instability does not start at a — 0 (Ai = 4e 2 ). This is because while energy from 
the ((f), A) sector does not change, the perturbations in A will affect derivative couplings 
between A and y, giving a positive contribution to the energy. Nevertheless, for any a > 0 
the strength of the instability is finite for A 2 = 7 = 0 , and so by continuity there will be 
some range of A 2 and 7 near zero for which the symmetric vortex is unstable. 

We begin by describing the perturbations. The cylindrically symmetric v<p = n solu¬ 
tion will have two independent functions which we call $ and v [[|: 

4>(r,6 ) = $(r)exp(in6 ) (2.3) 
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(2.4) 


^ n t \ 

A e = -v(r) 
r 


where r and 9 are polar coordinates on the 2-plane containing the vortex; we will consider 
n — 2 . As usual we have Aq — A 3 = A r — 0. We parametrize a perturbation about this 
solution as: 


(f>(r, 6) — Q(r)exp(2i6) + rj(f)(r, 9) 

(2.5) 

2 

Ae = ~(v(r) + rjv(r, 9)). 
r 

( 2 . 6 ) 

A r — rjw(r, 9). 

(2.7) 

where 7 << 1 is a small parameter. 

In particular, when a — 0, perturbations of the form 


(f>(r, 9) — (n$(r) 

( 2 . 8 ) 

Ai = + rjeijdjU). 

(2.9) 

where i, j — 1,2 and A^ is the cylindrically symmetric solution of the equation of motion, 
produce a zero mode if to satisfies Au — 24 )2 co. These zero modes of the energy may 
be resolved into eigenfunctions of the azimuthal angular momentum operator, so we may 
label the solutions as u: — u) m (r)cos(m9). For m — 1 the perturbation is nothing but a 
translation of the vortex in its own plane, while modes with m > 1 correspond to distortions 

of the string core. For finiteness of energy one must also require that m < where is 


the winding number of (f>. 

For a > 0, one can generalize these perturbations so that they reduce to the above 
form when a — 0. The translational zero mode exists in both cases; however, modes with 
m > 1 now decrease the energy, signalling the instability of the symmetric solution. 

We now consider this analysis for the combined ((f), y) system. Our boundary con¬ 
ditions require the winding numbers of cf> and y to be u w = 2 and u x — 1 respectively. 
The ((f). A) sector is identical to the system in [ 16 ] with — 2. It is sufficient to consider 
the m — 2 perturbation in the ((f>, A ) sector without disturbing the configuration of the 
scalar held y, since this only underestimates the extent of instability. Let us first expand 
the change in the vortex energy due to this perturbation in a power series in a and 7 for 
fixed A 2 , keeping k — 0 for now. We separate the change in energy into 7 -dependent and 
7 -independent parts: 

8E = 8E\ + 5E 2 (2.10) 
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where the first term on the R.H.S is 7 -independent. In fact, 8E 1 is equal to 8E for the 
system in | 6 ]. SE\ vanishes for a- = 0. so we can expand it in a power series in a; 


5E l = a Cl {K 0 + K t a + ....) ( 2 . 11 ) 

where the coefficients are all independent of the parameters of the x sector. 8E 2 is the 
7 -dependent part; this part is regular (at zero) in both 7 and a (regularity in 7 follows 
from the fact that 8E 2 is zero for 7 = 0 ). Therefore it can be expanded as 

5E 2 = 7 C 2 (Lo(«) + £i(«)7 + •••) (2-12) 

All the a-dependent coefficients of the expansion are regular in a. 

We have already seen that 8E 1 < 0 for a > 0 (and hence K {) < 0), while earlier 
we argued that for 7 > 0, 8E 2 > 0. There is then a curve in the (7, a) plane, passing 
through 7 = a — 0 and extending into the positive quadrant, along which there is a zero 
mode, i.e., 8E\ + 8E 2 — 0. This curve separates the positive quadrant in the 7, a plane 
into two parts with opposite signs for 8E , and represents the line of marginal stability for 
the cylindrically symmetric solution. Since 8E is negative in the part containing the axis 
7 = 0 , for every a > 0 there is a range of values of 7 for which any cylindrically symmetric 
solution to the equations of motion is unstable. However, it is clear that because of the 
field x the v,), — 2 string can not split completely into two well-separated strings, since 
isolated v ( p = 1 strings do not have finite energy per unit length. The ground state in this 
region of parameter space then must be a localized but asymmetric cosmic string. 

This proves the existence of binary strings in this model when k — 0. The restriction 
to k — 0 is not necessary, since just as we found a range of binary-string solutions in the 
a — 7 plane, we can continue this region of solutions into the third parameter-dimension, 
corresponding to the parameter k. Since the contribution of this term is nonsingular, 
for small values of k the cylindrical solution is still unstable. There is then no fine- 
tuning necessary to obtain this class of solutions, provided only that the desired pattern 
of symmetry breaking takes place. In particular, if we assume that the second symmetry 
breaking takes place at a lower energy scale, one does not expect the contributions from 
the lighter scalar to affect the question of whether the = 2 strings are stable, and so we 
would expect the binary string to be favored over a wide range of parameters, so long as 
a > 0 . 
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(Note that we could have added to the Lagrangian a term (4>*x 2 + c.c.) which removes 
the global symmetry allowing independent rotations of 4> and y; it can be similarly dealt 
with, and will be discussed further at the end of this section.) 

Our second example of binary strings confined by a long-range force is provided by a 
broken U(l) x U( 1) theory. This time we consider three charged scalar fields 0. x and r 
with charges ( 1 , 0 ), ( 0 , 1 ) and (A, 1 ) respectively. 

The lagrangian is 


L = 1 Wl 2 - \f^- + f\D\\ 2 - 1 + \\D0 - V. (2,13) 

F and G are the kinetic terms for the two gauge fields T /( and B /( . We also have, 


V = Ai 


vl) 2 + A 2 (|y| 2 - u 2 ) 2 + A 3 (|t | 2 - u 2 ) 2 + a(\^ - u 2 )(|y | 2 - u 2 )+ (2.14) 


, 2\2 


, 2\2 


- vim 2 - vi)+ k 2 (\ x \ 2 



We parametrize the ratio of the vev's by 


a2 _2_. 2 = ^3 

S .. 7.1 ' ..2 ‘ 


(2.15) 


Once again we first restrict the parameter space by putting = « 2 = 0. The vev’s 
of </>, y and r are respectively v i, 0- ! l an d 7 't’i- We expect to Hnd stable solutions only for 
boundary conditions for which the phase of each scalar field has an integer winding number 
at large distances; a minimal choice gives the winding numbers 1 , -1 and 0 , respectively. 
With this choice of couplings and boundary conditions, r will couple only weakly to the 
vortex, so we expect that |r| ~ everywhere. 

Let us momentarily consider the limit a = 7 = 0; the r sector contributes nothing 
to the vacuum energy, and the 4> and x sectors decouple. Two independent strings are 
formed, one coming from the winding of <f> an d the other from the winding of X- Although 
these two vortices will each be cylindrically symmetric and stable, the two-vortex state 
has a trivial zero-mode describing the relative displacement of the two vortices. This zero 
mode becomes an instability when one allows a to be greater than zero while keeping 
7 = 0, since now the vortex cores repel one another (with a short-range force). On the 
other hand, allowing 7 to be non zero while keeping a — 0 will also eliminate the zero 
mode, as the half-integer charge of r confines the singly-charged vortices in pairs with 
Qa + Qb = 0 mod 2 . For a = 0 the ground state will be cylindrically symmetric, but 
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taking a finite and increasing 7 from zero will clearly lead to an asymmetric ground state, 
since in that case the confining force is long-range but small. One can now argue (as in the 
first example of this section) that allowing both a and 7 to be greater than zero one will 
obtain a curve of marginal stability in the positive quadrant of the a, 7 plane; the part of 
this quadrant between the axis 7 = 0 and the curve of marginal stability is the region of 
parameter space for which doubly-charged, cylindrically-symmetric vortices are unstable 
to formation of confined pairs of singly-charged vortices, i.e. binary strings. 

Our initial restriction on the other scalar self-couplings is once again inessential. One 
can now add two more dimensions (corresponding to K\ and n- 2 ) to our parameter space, 
and the region of parameter space for which the binary string will be the ground state will 
be determined by some inequalities among the different contributions to the kinetic and 
potential energies. Since even a modest hierarchy among the rev's v t is sufficient to ensure 
that these inequalities pose little restraint on the parameter space, we see that for these 
models the binary string is technically natural and even generic, in the same sense that 
type I and type II vortices both occur in the phenomenology of superconductors. 

2.1. Goldstone modes and the binding force 

An interesting feature of the above models is the presence of a Goldstone boson in 
the broken phase. In both models we have omitted a scalar self-coupling which links the 
phase rotations of the different fields, and as a result there is an extra U(l) symmetry 
which is not gauged and is spontaneously broken along with the gauged U(l)’s. This extra 
global symmetry will be preserved by renormalization if originally present (unless, e.g ., 
there are fermions and the global symmetry is anomalous), and could arise as a result 
of discrete symmetries in a more fundamental theory, or through the absence of suitable 
renormalizable couplings in a unified model. In realistic models there are no massless 
scalars, and so either the symmetry-breaking coupling is present, or else perhaps the extra 
U(l) can be gauged. When the symmetry-breaking coupling is small enough the binary 
strings will exist by our earlier continuity argument. 

Let us define what we mean by the binding force. Consider quasi-statically separating 
the components of the string transversely along their length. Although these strings cannot 
be divided into their constituent strings, there will be deformations that stretch the binary 
string into two string-like boundaries of a domain-wall-like structure. Thus, consider a 
one-parameter family of deformations of the binary vortex in the first model, where the 
field © is constrained to vanish at two different points in the transverse plane at a distance 
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R. Finding the minimum energy solution to the equations of motion with this constraint, 
we define the binding force as (-) the derivative of the energy with respect to the parameter 
R. 

The fields forming the domain wall are massive in the true vacuum and the domain 
wall has a thickness h of order -A-. The energy of the domain wall is then a linear function 
of its area when the area is large compared to the thickness. For binary strings this 
translates to a binding force that is constant when R is large compared to the thickness 
of the domain wall. When the ‘domain wall’ is relatively thick {h » R) the energy 
grows less than linearly with the separation R. perhaps logarithmically. In fact, for small 
R. the fields in the region between the strings generally do not resemble a domain-wall. 
This small-separation case is quite typical of grand-unified models of binary strings, since 
the binding force and the short-distance repulsion may come from vastly different energy 
scales. Moreover, the nature of the binding force affects the action describing the motion 
of the binary string||. 

If the Goldstone mode is present, again the restriction to values of R much larger than 
h is insufficient. Nevertheless, by its nature the Goldstone mode is derivatively coupled and 
for large enough R the contribution from the massive fields to the domain wall energy will 
dominate the effects of the Goldstone boson due to their linear increase with R. Because 
of the strong changes in the fields near the string cores, the Goldstone boson may couple 
strongly to the motion and internal excitation of the strings. 

3. Molecular binaries 

Finally we consider models in which the attractive force between the two vortices 
or strings is short-ranged, and the existence of binary strings is due to a minimum in 
the vortex-vortex potential. The energy required to separate two vortices is then finite, 
but for strings of cosmological length, separation is energetically impossible. Consider a 
17(1) x 17(1) gauge theory with complex scalar fields (P\-02 and a neutral scalar a, with 
charges (1,0), (0,1) and (0,0) respectively. The Lagrangian is, 

L = \\D4>0 - + i|D^ 2 | 2 - i+ i|a<r| 2 - (3.1) 

where F and G are the held strengths of A and B respectively, D tl <f>i = (<0 M + ieA^cf) i, 
and = (dfj, + ieB^fa. The potential is 

V(0i,</> 2 ,cr) = Ai(|0i| 2 - juf) 2 + A 2 ct 4 + m 2 cr 2 - A 3 (p 2 - |</>i| 2 )a 2 + (3.2) 
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a 4 (|^| 2 - ti) 2 - X 6 (j 4 - \<h\ 2 )0*i - Wi\ 2 ) + Xe04 - l&| V- 

Note that the off-diagonal quartic couplings A 5 and Ag have the ‘wrong sign’; for a stable 
vacuum to exist this potential must be bounded from below. If the masses and coupling 
constants are all positive, a sufficiency condition for the positivity of the energy is easy to 
write: defining V — V' + m 2 a 2 we see that the positivity of V' is sufficient to guarantee 
the positivity of V. 

Now define x\ — /if — X 2 — cr 2 , X 3 = fi 2 — (fy, and in a matrix notation we can 
write V/ = X T AX ; where X = (X1X2X3) and 

/ Ai — A 3/2 —As/2\ 

A = I —A 3/2 A 2 Ag /2 I . (3-3) 

\ —A 5/2 Ag /2 A 4 J 

Vt is positive definite if A has no negative eigenvalues, which implies that the coeffi¬ 
cients of the characteristic equation of A must all be non-negative. This gives rise to some 
inequalities involving the couplings: 

TvA — Ai T A 2 T A 4 >= 0, (3-4) 

detA = 4A1A2A4 + A3A5A6 — AiAg — A2A5 — A3A4 >= 0, (3-5) 

4[AiA2 + A2A4 + A4A1] — A3 — A5 — Ag >= 0. (3-6) 

When these conditions are satisfied the potential energy can be written as a sum of 
squares by diagonalizing A so it is explicitly non-negative, furthermore, we have arranged 
that V — 0 in the desired vacuum. We will assume that all the A’s are positive thus 
satisfying (|T4|)at once. We shall come back to the other conditions on the A’s later, when 
the other requirements are known. 

First let us note that with this Higgs potential, the U(l) x U{ 1) symmetry is sponta¬ 
neously broken; the vacuum expectation value of the symmetry breaking Higgs fields are 
\(f)\,vev \ — l J i and \ 4 > 2 ,vev\ = P' 2 ; and the vev of the singlet field a is zero in the true vacuum. 
The manifold of vacua has nontrivial topology; 7 Ti(f 7 (l) x 1/(1)) = Z x Z implies the ex¬ 
istence of two kinds of strings: those with a flux of A and those with a flux of B. We shall 
call them type A and type B vortices, respectively. These can occur independently , so the 
binding which occurs is through short-range forces only. As we shall see, the vortex-vortex 
binding can occur between vortices of different types. 
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Vortices of circulation one in either A or B alone are always stable. We are interested 
in the case of a pair of parallel A and B type strings with a vev of a in the core of the A 
type string. Let us review the conditions for which a would get a vev in the core of the A 
type string. The terms in the potential involving a are, 

V(a) = X 2 cr 4 + m 2 a 2 - A 3 (p 2 - 4>\)(r 2 . (3.7) 

At the core of A type string we can write a local potential for a with other fields replaced 
by their values at r — 0 , (j>\ — 0 , (j )2 — ^ 2 - 

V(a) — X 2 cr 4 + m 2 a 2 — Xs(^i)cr 2 . (3-8) 

There is a non trivial minimum of this potential at a 2 ev — m , when A 3 P 2 > m 2 . 

We can write a similar potential at the core of a string of type B; however, the sign 
of the term with A 6 in the potential has been chosen so that the minimum of the ‘core’ 
potential is now at a vev — 0. The true configuration of the a held inside the A vortex is 
of course dictated by a balance between the gradient and potential contributions to the 
energy. However, when the condition 


\ 3 l4 > m 2 


(3.9) 


is fulfilled, a nonzero value of a vev is preferred at the core of the A type strings because of 
the gain in potential energy. 

The conditions for obtaining a condensate of a ‘spectator’ scalar held were examined 
in let us briefly recapitulate. We begin with a static, cylindrically-symmetric vortex 
solution (Nielsen-Olesen vortex) for ©\ and the gauge held A /( . neglecting ©2 for now. We 
wish to see if a — 0 is indeed the vacuum state; we examine small oscillations of a around 
the background a = 0, a = ao(r)e lujt . The linearized equation of motion of a is, 


dV 0 (r) 

dr 2 


1 da 0 (r ) 
r dr 


Pa 0 (r) = u 2 o 0 (r) 


(3.10) 


When there are solutions of ( |3.10|) with tv 2 < 0, a — 0 is an unstable solution of the clas¬ 
sical equations 0. When m 2 = 0, the potential energy of this 2d Schroedinger equation, 
P = A 3 (p 2 — d>\) is everywhere negative when the background is the Nielsen-Olesen vortex 
solution (f)i(N.O). The argument in [| 2 |then proceeds that, since a negative definite po¬ 
tential in two dimensions always has a bound state, uj 2 has a negative value for m 2 = 0 ; 


fO 







by continuity, then, in the neighborhood of m 2 = 0 , a — 0 is an unstable solution which 
relaxes to a non zero vev for a. Indeed, a develops a nonzero vev whenever m 2 < /i 2 , where 
—p 2 is the most negative eigenvalue for to 2 in equation ( |3.6| ) with m — 0. 

Although we cannot express /i directly in terms of the parameters of the theory, it 
can be easily approximated by doing a simple variational calculation, or even more simply, 
estimated by making an ansatz for the zero mode, parametrizing the shape of the vortex 
by elementary functions. Choosing the latter method (taking the ansatz for ay to be given 
by the difference of two exponentials) we get 

~ 2.4/c 2 -0.6A 3 /zf. 

where k = \/‘2Xi // 1 . The first term on the RIIS is from the kinetic energy and the second 
is from the potential energy of the a condensate. There are two things to note: (a) This 
value of cu 2 underestimates the actual value of // 2 . However, we can also obtain an upper 
bound; for a vortex with |</>i(r)| that monotonically increases from (pi — 0 at r = 0 to 
(pi — Hi at r = 00, we must have ji 2 < A 3 /U 2 . (b) The value of p 2 is an increasing function 
of the range of the potential, when the strength of the potential is held constant. This 
is because the potential energy term P(r) increases in magnitude point-wise as the range 
of the potential is increased. In our case the range of the potential is 1/k, so /i 2 is an 
increasing function of k. 

We want to show that there is a range of parameters for which (i) a develops a vev 
in the core of the A type strings; and (ii) the asymptotic forms of the (pi and a fields is 
\(pi | (r) ~ iii — (pie~ kr and |er|(r) ~ aoe -mr . (i) will be true if the condition m 2 < p 2 is 
satisfied. Since /i 2 < A 3// 2 we can define a vortex parameter /y, (depending on the vortex 
field background (pi ) by 

Im 2 Ui =(3.11) 

with 0 < <= 1 . 

Condition (ii) for |<r|(r) simply follows from the linearized classical equations in the 
asymptotic region. Similarly we obtain the correct asymptotic behavior for <pi if m 2 > 
8 Ai/i 2 and e 2 > 2A 2 . For consistency with (i) we then must have, 

f M a 1 )A 3 >8A 1 . (3.12) 

For given values of f.i 2 and A 3 there exist a range of values for Ai for which this consistency 
condition is satisfied. This is quite natural since the width of the Nielsen-Olesen vortices 
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increases as Ai —> 0; in that limit, ft 2 tends to its maximum value: /r —► A 3 // 2 . Thus, given 
any A3 > 0, there is a range of values of Ai extending to zero satisfying /^(at.o)^3 > 8A1. 

Having established the asymptotic behavior of the (f> 1 and a fields when A 3 > Ai we 

now examine the behavior of a B-type vortex in the background of an A-type vortex, which 

as we see has a a condensate in the core. The scale that determines the width of the B type 

vortex is the mass of the (f >2 field in the true vacuum, width(B) ~ , 1 . We assume that 

2M2 

the ((f) 1 , a) and the (f) 2 sectors have a large difference in mass scales, /r| >> /if ~ m 2 . This 
would be natural for A 5 = A 6 = 0, when the two sectors are decoupled. In a GUT version 
this hierarchy would result from symmetry breaking at different scales; in the present case, 
it does not need to be very large. In this simplified model we assume A 5 , A 6 to be relatively 
small: A 2 , A 3 , A 4 >> A 5 , A 6 - 

The large mass scale of the B type vortex ensures that the B type strings are very 
thin compared to the A type ones. We will then approximate the ‘shape’ of the B vortex 
effectively by a delta function, with the field (f >2 dropping abruptly to zero from its true vev 
H 2 at the core of the vortex, for the purpose of considering the interaction of the vortices. 
If we bring such a (thin) B-type string near an A-type string, while keeping them parallel, 
the structure of the vortices will not be changed appreciably by their mutual interaction. 
Indeed the interaction is purely through quartic Higgs couplings of A 5 and Ag, which can 
be treated as a perturbation. 

We can now calculate the interaction energy per unit length of a pair of parallel A 
and B type strings, V in t = ^l' 2 [X^a 2 (r') — A 5 (/i 2 — 4>\(r'))\, where // ~ / 1 2 is determined by 
integrating the interaction energy over the cross section of the thin B-type vortex, and r' 
is the position of the B-string relative to the core of the A-string. 

When the B string is far away from the core of the A string we can use the asymptotic 
forms of the (f> 1 and a fields: 


\(f>i(r)\ ^ Hi - (pie v 8AlMl(r) ; 

|cr(r)| ~ ae ~ mr ; 

V int = p 2 [A 6 a 2 e- 2mr - 2A 5 /xi0e-^ Al ^)] + h.o. 

At r —> 00 , dV Q^ is positive (m 2 > 2Aipf), giving an attractive force between the vortices. 
This force decreases to zero and changes sign at the equilibrium point r = r eq which is 
given by 2X 6 a 2 me~ 2ujre ‘i — 2 A 5 pi</>ie - v /AlMl(re ‘ i ) + h.o. = 0 . 
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The vortex-vortex interaction leads to a equilibrium separation of the strings, r eq . By 
choosing an appropriate ratio of the constants A5 and Xq, r eq can be made to fall within 
the range of validity of our approximations. In the present case, A 6 must be much larger 
than A5. (That is, we have balanced a somewhat-longer range attraction with a stronger, 
but shorter ranged repulsion.) This will obviously lead to a binding energy per unit length 
between the two strings that is several orders of magnitude smaller than the mass per 
unit length of the binary string, since the attractive force is proportional to A5. To get a 
larger binding energy one could increase A5 (and then A (i as well), in which case one must 
replace the asymptotic forms of the (f> 1 and a fields by the exact classical solution, since 
r eq moves into the core of the A-string. There is a priori no reason that the qualitative 
features would change drastically for more strongly bound strings; there will be some range 
of couplings for which the binary-string ground state is obtained, with the limit on the 
strong-binding side being set by the fact that as the couplings get stronger and the strings 
are more closely bound, it becomes energetically favorable to deform the A-string and 
obtain a cylindrically-symmetric configuration. In any case we have an existence proof 
of vortex-vortex binding at finite separation in this model, and presumably others of the 
same class can be found. 


(Returning to the constraints on our scalar potential, the nontrivial conditions ( 3 . 5 ), 
|) and (| 3 . 12|) in the above model are met for a wide range of (positive) parameters A. 
For instance when A5 and A 6 are sufficiently small the conditions essentially reduce to 

4A1A2A4 — A3A4 > 5 , 

4 [AiA 2 + A2A4 -f- A4A1] — A3 > 5 , 


/<£A 3 > 8A1, 


where 5 is a small positive parameter. If Ai = efy A3, with e < 1 , then conditions (| 3 . 5 |), 
Q 3 .G|)and ( ]). 12 | ) are easily satisfied by any partial ordering A2,A4 > A3 satisfying 4A2 > 

A3 /eU-) 


To compare this with the previously-constructed models, the most significant differ¬ 
ence is the finite binding-energy per unit length, which is much smaller than the string 
tension. The individual strings in this binary are independently stable; if they were initially 
separate, they would be very unlikely to combine because the interaction energy between 
two intersecting segments of string could not compete with the kinetic energy of the rest 
of the string under normal circumstances. On the other hand, if such a binary string were 
formed it would be very difficult to dissociate, since a macroscopic string would have a 
very large binding energy, and it would be necessary to transfer sufficient energy to the 
whole string, rather than an individual segment. 
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4. GUT considerations 


Binary strings of the ‘confining’ type will occur in grand-unified theories provided two 
conditions are satisfied: the sequence of symmetry-breaking phase transitions gives rise to 
cosmic strings and broken discrete symmetries in the necessary order, and the coupling 
constants are such that multiple-winding strings are unstable. The first of these conditions 
is the more stringent one. The question of whether the pattern of symmetry breaking is 
appropriate depends on the representations of the symmetry-breaking order parameters 
(e. g .Higgs fields) and the scales at which the symmetries are broken. (Note that although 
we will be describing the symmetry breaking in terms of elementary scalar fields, there 
is no reason why other symmetry-breaking mechanisms, such as fermion condensates or 
Wilson lines, could not play a role). 

The first type of confining binary, in which two similar strings are bound together 
by a (vestigial) domain wall, is the most generic, since it could occur in any grand-unified 
model in which strings are first formed and at a later phase transition bound domain walls. 
It may be surprising that a stable remnant persists in this case; however, we observe that 
the orientation of the domain wall is not gauge-invariant, but changes direction under 
the Z 2 discrete gauge symmetry spontaneously broken by the lighter scalar field, so the 
string can assume either orientation as the boundary of the domain wall. (In terms of 
the holonomy of the gauge field, the same group element, (- 1 ) in the representation of the 
charge -1 field, is obtained in passing around the string in either direction, and likewise 
passing in either direction through the domain wall.) In the next section, we remark on 
how this circumstance allows (n= 2 ) strings, whether binary or not, to survive the phase 
transition in which the domain walls form. 

This type of confining behavior can occur in a well-known model of SO(IO) grand- 
unification described by 0 in their analysis of string-bounded domain walls. In that 
symmetry breaking sequence, 

Spin(10)-^Spin(6) 0 Spin(4)™SU(3) 0 SU( 2 ) 0 U(1)^*SU(3) 0 U( 1) (4.1) 

cosmic strings are formed (of the ‘Alice string’ variety |||] H) at the first phase transition; 
the 126 vev is not C-invariant, and so at the second phase transition the strings bound 
domain walls. The surviving strings will be of the binary or the n — 2 type, depending 
on the relative sizes of the scalar and gauge couplings, as previously discussed. Other 
examples may be similarly constructed. 
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Incidently, this same model (in a different region of parameter space) might exhibit- 
the ‘molecular’ variety of binary string. At the phase transition in which the domain 
walls are formed, there are independently formed [7|J stable cosmic strings (corresponding 
to the discrete gauge symmetry ( —l) Fl ° G Spin( 10)). They might be expected to be 
comparable in abundance to the remnant n — 2 strings remaining after the domain walls 
have disappeared. One could then extend the scalar sector of this model so that the 
situation of section 4 is realized, that is, a weak attraction at long distances between the 
strings, and a short-distance repulsion, due to a scalar condensate in the core. In view of 
the somewhat dubious cosmology of this type of defect, we will not develop such a model 
in detail. 

If we consider grand-unification based on the Eq gauge group, there are more possi¬ 
bilities for the unbroken discrete gauge symmetries and hence a variety of string solutions. 
So, for instance, if at the highest scale we break E G to SU(5) 0 1/(1) 0 U(l) via one or 
more adjoints (78's), we can obtain independent cosmic strings at subsequent stages of 
symmetry breaking, either in combination with or independent of SU( 5) breaking. These 
could furnish examples of any of the varieties of binary string we have described. 

We will present one concrete example along a somewhat different line. Consider an 
Eq GUT model, with a symmetry-breaking sequence 

E g ^SO{ 10) 0 Z^SU{ 5) x Z 2 0 Z 2 ^SU{ 5) 0 Z 2 -> ... (4.2) 

(we leave the subsequent symmetry breaking to the reader’s imagination.) The represen¬ 
tations of the scalar vev’s are, respectively, a 27 0, a 351' y, and another 27 r. (Note that 
the 351' is in the symmetric product of two 27’s, and the direction of the vev is chosen to 
lie in the 126 direction in its S'O(IO) decomposition, so the Z 2 string is similar to those in 
S'O(IO) GUT models). Our nicknames for the scalars is meant to suggest the similarity to 
the second model of section 2; indeed, the only difference is that in this case the discrete 
symmetry group has a different breaking pattern, Z 4 0 U (1) —> Z 2 0 Z 2 —► Z 2 . Because of 
the Z] breaking the composite strings may form at the second phase transition, as well as 
at the third. 

Finally, we once again emphasize that the point of this exercise is not to find compelling 
grand-unified models with binary strings, although the first example would do for that, 
but rather to show that in non-minimal unification the conditions for the existence of 
composite strings are easily satisfied. 
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5. Formation and Evolution of Binary Cosmic Strings 


The formation of binary cosmic strings begins by the Kibble mechanism [|10]|: the 
Higgs scalar develops an expectation value below some critical temperature T c , which is 
disordered on scales larger than some correlation length determined by the kinetics of 
the phase transition. As the scalar field relaxes, defects are formed at the boundaries of 
domains with different orientations of the condensate. At lower temperatures, the scalar 
field fluctuations are suppressed and the network of defects begins to evolve. 

The details of binary string formation depend on the mechanism responsible for the 
internal structure and the evolution subsequent to the initial formation of the strings. In 
the case of confined binary strings domain walls are formed at a later phase transitions 
when one of the discrete gauge symmetries is broken. Since the strings arising from the 
earlier symmetry breaking are charged under the discrete symmetry, they will become 
boundaries of domain walls. As the network of domain walls forms in the presence of the 
strings the relative orientations of different components of the boundary of a given wall is 
determined by chance. Then, as the network evolves, some components of the boundary 
will meet as the domain wall shrinks; those that have opposite orientations will annihilate, 
while those segments of string that have the same orientation will form linear defects with 
winding number n = 2. These will be binary strings, when the couplings are such that 
the cores are repulsive; otherwise they will form metastable, cylindrically-symmetric n — 2 
strings. (In both cases the strings are unstable to decay by Z 2 monopole-pair production 
ITU, through a greatly suppressed tunneling process; because of the separation of the cores 


the rate will be even more heavily suppressed in the case of binary strings.) 

The situation is similar with the other confining model, in which two different types 
of string are bound together at the discrete-symmetry-breaking phase transition. Here, 
however, the domain-wall can end with either orientation, on either of the types of string. 
Generally we may assume one of these strings is much lighter and more abundant than 
the heavier ones initially formed. It is then reasonable to expect that while many of the 
lighter strings may annihilate, nearly all of the heavier strings should be paired up with 
the more abundant lighter strings. (Note that in this case if the couplings do not make 
the cores mutually repulsive, the final state will have strings with a gauge-field flux being 
a composite of that comprising the initial strings.) 

The circumstances under which the ‘van der Waals’ or ‘molecular’ strings may be 
formed cosmologically are quite unclear. There may be circumstances under which pre¬ 
existing strings from an earlier phase transition could serve as nucleation sites for the 
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subsequent transition or otherwise influence the ordering of the scalar field condensate, 
in such a fashion that the binary strings are formed together. Otherwise, since the sepa¬ 
rate strings are stable and the interaction in this case is short-ranged, there is vanishing 
probability of forming binary strings when the component strings are formed separately. 

After formation the evolution of a network of binary cosmic strings is essentially the 
same as for ordinary cosmic strings, since at large distances the string mechanics is domi¬ 
nated by the string tension (which is nearly the sum of the tension for each component). 
One difference is that the initial network of strings formed by the collapse of the domain 
walls may be somewhat different from the usual string network, but differences of this 
type will presumably be washed out by the subsequent evolution of the string network. 
Another effect (albeit a small one) which should be noted is that the effective string tension 
and mass on cosmological scales |l2j| will be affected by the internal degrees of freedom 
(twisting modes). 

It is more interesting consider the coupling of the internal modes of the binary cosmic 
string to the various epiphenomena of cosmic string solutions: superconductivity, fermion 
zero modes, the kinks and cusps of classical string evolution. Some of these questions will 
be addressed in the companion paper [0. 
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